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Abstract 



^ I Low energy neutrinoless lepton flavor violating (LEV) processes are studied 

• I in an extension of the Standard Model (SM) by heavy SU{2) x U{1) singlet 

PsJ ' Dirac neutrinos. An upper bound procedure is elaborated for the evaluation 

^— ^ ■ of amplitudes. A comment on the extraction of heavy neutrino mixings from 

Q^ . astrophysical observations is given. For processes not treated in the applied 

^^ \ model the formalism for evaluating the branching ratios (BRs) is presented. 

The processes previously studied in the model are carefully examined. Some 
of the previous results are improved. Special attention is paid to the structure 
of the amplitudes and BRs as well as to the relations between BRs of different 
^^ . LFV processes. Numerical analysis of the BRs is done. The decoupling of 

heavy neutrinos is discussed and it is explicitely shown that the very heavy 
neutrinos decouple when the upper bound procedure is applied. The upper 



;h ' limits of the BRs are compared with the current experimental upper bounds 

and the processes interesting for the search for LFV are proposed. The LFV 
decays are shown to be unsuitable for finding upper bounds on "diagonal" 
LFV parameters. The 5-meson LFV processes are suggested for the search 
of LFV in future S-factories. 
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I. INTRODUCTION 

If the instanton effects |[I[] are neglected, lepton flavor and lepton number are conserved 
in the Standard Model (SM). Recently found atmospheric neutrino oscillations indicate 
that neutrino masses are nondegenerate and the lepton flavor is not conserved. Independent 
confirmation of the deviation from SM is expected to manifest as a nonconservation of lepton 
flavor/number (LFV/LNV), as a breaking of lepton universality, in CP violating processes 
which are not consistent with SM etc. 

The problem of LFV/LNV is related to the physics beyond SM and includes various 
physical areas [^ : atomic physics (e.g. muonium-antimuonium conversion), nuclear physics 
(yU — > e conversion, double beta decay), low energy hadron physics (leptonic and semileptonic 
decays of mesons and tau lepton), problem of CP violation etc. 

LFV have been found in various extensions of SM 0-^. Here, LFV is studied within one 
of two extensions of SM by heavy neutrinos 0J^ , obtained by adding additional heavy Dirac 
neutrinos to it. It is referred here as the V model [§]. Due to the Dirac character of the 
heavy neutrinos, there are no LNV processes in this model. The other model 0, obtained by 
extending SM with additional heavy Majorana neutrinos, has some renormalization problems 
and light neutrino mass problems 0. Besides the additional heavy Dirac neutrinos, the V 
model contains three massless neutrinos. It should be noted that in this work the V model is 
used phenomenologically. Any model with the same gauge properties and about equally large 
heavy neutrino masses would give the same results, regardless whether the light neutrinos 
are massless or have masses which are in accord with the present experimental data. 

The extensions of SM by heavy neutrinos contain a Cabibbo-Kobayashi-Maskawa type 
matrix for leptons (LCKM). In general, the elements of this matrix are not known. Experi- 
mental and theoretical constraints exist only for some specific sums of the matrix elements of 
the heavy neutrino part of the matrix. Therefore, the LFV amplitudes cannot be evaluated 
exactly, but only the upper bounds on their values may be found [0. The evaluation is 
especially complicated when the amplitudes contain expressions with more than two LCKM 
matrix elements. In this paper a method for evaluation of the upper bounds of amplitudes 



found in the previous publication [|TU| is improved. The method gives upper bounds for 
all values of the model parameters, but in some directions of the parameter space it is 
not very restrictive. It is explicitly shown that the upper bound procedure leads to the 
decoupling of the heavy neutrinos in the infinite mass limit, showing that the "nondecou- 



pling" of heavy neutrinos ||TT| , [12[| is only a transient effect, appearing with enlargement of 
the heavy neutrino mass. It should be noted that this "proof of generalization of the 
Appelquist-Carrazone theorem is based only on the requirement that the physical system 
can be described pertubatively, and is independent of the introduction of somewhat unde- 



termined maximal S'f/(2)2,-doublet mass term as in Ref. []12|. To give the feeling how large 
error can be introduced using the upper bound procedure elaborated here, a few branching 
ratios (BRs) obtained by the upper bound procedure are compared with BRs obtained using 
"realistic" LCKM matrices. 

The LFV processes are not very usefull for deriving upper bounds on the matrix elements 
of LCKM matrix. The amplitudes for these processes are proportional to the sums of 
products of the LCKM matrix elements and functions of heavy neutrino masses. Using 
the freedom to choose unknown phases of the LCKM matrix and heavy neutrino masses. 



these sums can always be set to be equal zero, even if the absolute values of nondiagonal 

elements of the LCKM matrix are different from zero. The present limits on the LCKM 

matrix elements are derived from the measurements of lepton flavor conserving processes 

13| , more precisely, from the estimates of deviations of the corresponding decay rates from 



the SM results. For each row (a row corresponds to a specific lepton /) of the LCKM matrix, 
these data give a limit on the sum of squares of absolute values of the matrix elements 
corresponding to the heavy neutrinos, (s^)^. Knowing the upper bounds on (s^')^-s one 
may derive the upper bound for BR of any LFV process. One of the aims of this paper 
is to derive the upper bounds of BRs for all low energy LFV processes in the V model. 
The processes having comparable theoretical and experimental upper bounds of the BR, 
or theoretical upper bound larger than the experimental one, are interesting for further 
experimental investigation. 

Neutrino oscillations of two massless neutrinos in supernovae have been shown to give 
a very strong upper bound of two of the LCKM matrix elements in the part of the ma- 



trix corresponding to the massless neutrinos [14|. Here, the analysis has been repeated for 



three neutrinos, hoping that the upper bounds for other "massless neutrino" LCKM ma- 
trix elements may be derived. The knowledge of nondiagonal "massless neutrino" LCKM 
matrix elements may, in principle, lead to better upper bounds on some combinations of 
"heavy neutrino" LCKM matrix elements than the terrestrial experiments. Unfortunately, 
the analysis made here shows that the three-neutrino oscillations do not give new constraints 
on any combination of "heavy neutrino" LCKM matrix elements. The only new information 
it gives is that the mixing between massless "mu" and "tau" neutrinos is smaller than the 
value obtained from the analysis of Super-Kamiokande data [Q, in which "mu" and "tau" 
neutrions were assumed to have small masses. 

Untill now, many of the low energy neutrinoless LFV processes have been investigated. 
Some of them were examined only within a few models, for instance LFV decays of heavy 
mesons. The neutrinoless LFV decays of B-mesons were studied in the frame of SM with 
additional Higgs doublet |TH], while the neutrinoless LFV decays of D- mesons have been 



studied in the frame of leptoquark models [T^ and a flipped left-right symmetric model |]T7 



Here, they are analyzed in the V model. Some of the low energy LFV processes have not 
been studied in the frame of the V model. Among them are the muonium-antimuonium 
(M <-i> M) conversion and neutrinoless LFV violating decays of the Z boson. The results are 
also given here. Some of the neutrinoless LFV processes have been analyzed in the V model, 
but the analysis is incomplete |T^ or there are some errors in the expressions for amplitudes 



or decay rates [p!0|-[l^. Here only the corrections to the previous results are given. 

On the quark and lepton level there are only a few Feynman diagrams (composite loop 
functions) that contribute to any neutrinoless LFV decay amplitude. If two neutrinoless LFV 
processes contain only one common composite loop function, the ratio of corresponding BRs 
is independent of the \^-model parameters. Therefore, roughly speaking, knowing one BR, 
BRs of processes comprising the same basic Feynman diagram may be evaluated without 
the knowledge of parameters of the V model. If LFV decay amplitudes contain different 
loop functions or more loop functions, the ratio of the BRs depends on V^-model parameters. 
Nevertheless, the mass dependence of the ratio of the BRs simplifies in the limit of large 
heavy neutrino masses. Most of the amplitudes become dependent essentially only on one of 
the composite loop functions. In that limit, the ratios of the BRs having the same dominant 



composite loop function become independent of the V^-model parameters. Experimentally, 
for most neutrinoless LFV processes, only the large heavy neutrino mass limit is interesting, 
because, with few exceptions, only in that limit BRs assume the values comparable with the 
present day experimental limits. A comparitive analysis of the amplitudes and BRs of all 
neutrinoless LFV processes is presented. 

In the section || some properties of the V model, relevant for further discussion, are given. 
Discussion on the limits of the model parameters is given in the section |T|. The amplitudes 
for the neutrinoless LFV processes not studied in the V model, and some improvements 
and corrections of the previous results are presented in section |V]. The amplitudes and 
BRs of LFV processes are studied in section |V]. The numerical results and comparison with 
experimantal limits are also given in section 0. Conclusions are summarized in section ^1. 



II. COMMENTS ON THE MODEL 

Here, a model with additional SUl(2) x f/(l) singlet Dirac neutrinos, which have large 
mixings with the SM leptons, is used in the calculations. The masses of the singlet neutrinos 
are not restricted by the SUl{2) x f/(l) breaking scale. The large mixings and the large 
masses are necessary conditions for obtaining observable LFV decay rates. 

In the model considered here P,pD|-p2|, the total lepton number (L) is conserved. For 



each SM neutrino one left-handed and one right-handed singlet neutrino is added, although, 
in principle, the structure of the mass matrix permits addition of an arbitrary number 
of pairs (n^j) of left-handed and right-handed neutrinos {Vnji models). Lepton number 
conservation gives such a structure to the mass matrix which automatically leads to three 
massless neutrinos at any order of the perturbation theory ||19|| . 

Since the new neutrinos are SU{2)l x U{1) singlets, the structure of the lepton interaction 



vertices in the weak basis is the same as in SM |n]. However, in a transition to the mass 



basis, nondegeneracy of the neutrinos leads to the Cabibbo-Kobayashi-Maskawa (CKM) 
type matrix (Bin) in the charged current (CC) nlW vertices. As only a part of the mass- 
basis neutrinos interact with the Z boson, neutral current (NC) nnZ vertices {n is neutrino 
field in the mass basis) are also not flavor-diagonal, and contain matrix elements of the 
nondiagonal matrix (C„„). The NC IIZ vertices and the quark vertices are the same as in 
SM. 

The C matrix from the neutrino NC vertex may be expressed in terms of B matrices 
from the CC lepton vertex. Therefore, besides the SM parameters, the model depends only 
on the B matrix (or more precisely on the parameters defining the B matrix) and on heavy 
neutrino masses. The matrices B and C satisfy a set of relations stemming from the gauge 
structure (see e.g. |TT|), 

na+riR nc+nu 

k=l k=l 

nc+nR riQ 

/ , BikCki = Bii, ^Bij^Bij = Cij. (1) 

fc=i 1=1 

From the orthogonality relations for Bin matrix elements, phase arbitrariness of leptons 
and SU{nji) invariance of massless neutrinos lead to ncnn independent angles and [ug — 



1)(^_R — 1) independent phases of the B matrix [^,Q. Experimentally, only nc parameters 
s^' may be estimated. Therefore, the B matrix elements are undetermined even for the 
simplest case with two additional heavy neutrinos, nn = 2. Since the B matrix elements 
are unknown, the amplitudes of LFV processes cannot be evaluated exactly, but only upper 
bounds of the amplitudes may be found. One should mention that there exists a model 
for which amplitudes of LFV processes can be evaluated exactly, in the case oi ur = 2 
1^. Unfortunately, as mentioned before, it is excluded because of some renormalization and 
light-neutrino mass problems. 

The degeneracy of massless neutrinos allows one to write the B matrix in the following 
form [|,|19|,|5[ 



B,,^ = [{UDa)i,,^,{UGU,] 



k 



I), 



(2) 



where f/ is a unitary matrix. Da is a diagonal matrix and G is a matrix satisfying D'j^+GG'' = 
1. Indices i and I denote massless and massive neutrinos, respectively. From the structure 
of the B matrix, it follows that the massless neutrino CC in principle is not diagonal, leading 
to LFV [0,^,Q and nonortogonal effective weak-neutrino states ||26|, although neutrinos 
are massless. On the other side, the massless-neutrino NC, which contains the C matrix 
elements, is diagonal [|19|. Since there are no tree-level flavor violating neutral currents 
(FCNCs) in the massless neutrino sector, the universality of massless neutrino couplings is 
not satisfied, because, in general, the elements of the diagonal matrix Da are not equal. The 
nonuniversality of these couplings may have some astrophysical implications. 

As mentioned in Introduction, the B matrices are used to define the parameters s^, 
which are a measure of the deviation from SM, in the following way |ll2| , p7| -pO[] 



-.''l\2 



riR 



^BiN^B^ 



IN, 



(3) 



Because the definition of [s^Y contains Bi^^f matrix elements of the same lepton flavor, the 
term "diagonal" mixing(s) will be sometimes used in the text below. 

III. LIMITS ON THE MODEL PARAMETERS AND METHODS OF 
EVALUATION OF AMPLITUDES 

A. Experimental limits 



The parameters (s^)^ have been determined from the global analysis of the low energy 



tree level processes |]13| , p7| -p0[] . In these processes heavy neutrinos may manifest only in- 



directly, through change of the light (massless)-neutrino couplings. These couplings attain 



additional c2 factors, where {c^Y 



l-(s 



•^l\2 



Y^i=i BiuiB^i,.- The changes of the couplings 
could show up as a nonuniversality of CC couplings, as a deviation from unitarity of the 
CKM matrix, as a change of the invisible width of the Z boson etc. [0,Q. The best limits 
on the mixings s^'. 



[s^^f < 0.0071, 



[sl^f < 0.0014, 



< 0.033(0.01), 



(4) 



1^. The value in the brackets is vahd for SU{2)l x f/(l) singlet heavy 



were found in Ref. 
neutrinos. 

The more stringent limits on the Bi^ matrix elements were searched for investigating 
the loop effects in the lepton-conserving and lepton-violating processes. Direct limits on the 
parameters s'2 are not possible as the expressions derived from the loop amplitudes, which 
are constrained by experimental data, depend not only on the s'2 parameters but also on 
the Bin phases and masses of heavy neutrinos. Lepton conserving processes including heavy 
neutrinos in loops were studied by Kalyniak and Melo pl| , |32| . They studied the loop effects 
of heavy Dirac neutrinos on muon decay, universality-breaking ratio in Z —>■ II decays and 
Ar quantity. They found no new constraints on the s'2 parameters. The flavor nondiagonal 
(LFV) processes without light neutrinos in the final state were studied extensively both 
theoretically P,p|, p!T] , p!9| , pO| , p3| , p3| , p5[ and experimentally [0,^^-|3^. The advantage of these 
processes is that their observation would be a clear and unambiguous signal for LFV. These 
processes proceed only through loops. Using independence of the loop functions on the 
light neutrino masses and the orthogonality of rows of B matrix, the amplitudes of these 
processes may always be expressed in terms of heavy neutrino contributions only. Three 
of these processes, fi — > 67, fi —>■ 3e and e-fi conversion in Ti, gave new very stringent 
constraints on a specific combinations of heavy neutrino masses and matrix elements B^n 
and B,,N IIT2|. Particularly, near-independence of the fi —>■ ej amplitude on heavy neutrino 



■>fiN 



masses enables one to find the following very stringent mass independent limit. 



nR 



J2b;nA^^ < 2.4x10- 



(5) 



i=l 



No other constraints independent of heavy neutrino masses were derived from the LFV 
processes. It should be noted that the limit (j^) does not necessarily lead to new limits on 
the s2 parameters. The sum in (^ may be written in terms of the parameters s'2' and s'^ 
and a complex "cosine" of the "angle" between vectors {B^^-} and {B^Ni}, 



riR 

E 

i=l 



^liNiBeNi 






(6) 



where x°g 



X]"=i B*^NiBeNi/ s''2^ s'2; . Obviously a reduction of x°g may assure the fulfillment 
of the inequality (|^) without reducing the s2 parameters. Within the V model the explicit 



BfNjs^ 

ut r 
estimates of BRs for the processes including more than two Bi^ matrices were given for the 



first time in Ref. |[T0[| . 



B. A comment on astrophysical limits 



The masslessness of "light" neutrinos in the V model lead to the limits on some Biy matrix 
elements which can be derived from astrophysical observations. Valle and collaborators have 
noticed that the measurements of neutrino flux from the supernova SN87 leads to two very 



small lepton-massless-neutrino mixings ||14 



IS, 



\BruA < 10" 



(7) 



The result (|^) follows from an estimate of the v^-Vr conversion probability in the V model. 
To find whether similar upper bounds can be found for other massless-neutrino B matrix 
elements, their calculation is repeated here for three massless neutrinos. The motivation 
for such calculation is following. Through the orthogonality relations for B matrix elements 
(0), very stringent limits on the matrix elements Biy would lead to better upper bounds on 
nondiagonal mixings Y^=\ BiNiBpj^. than those obtained by terrestrial experiments. 

The derivation of the limits (^ is based on an analysis of neutrino oscillations of the two 
neutrinos for which the experimental upper bounds on s^ parameters are the weakest - s'^ 
and s'l . The oscillations of massless neutrinos are a consequence of an interplay between 
the CC and NC neutrino weak interactions |^ . They appear only if the universality of the 
NC interactions is not fullfilled and if the nondiagonal CC currents are different from zero. 



Following the notation of Refs. 
small parameters hi (for small hi, 



the deviation from the universality is described by 
hi ~ s2)- The massless-neutrino part of the B matrix is 
parametrized by one mixing angle 6, which is assumed to be small. The resonance condition 
reads 



2K 



/i? 



h"^ 

"'e 



l + hl 



(8) 



where Fg = ne/ijie + n„), F^ = 1 — Fg and rie and Un are the electron and the neutron 
number densities. As the experimental limits on he and h^- are much smaller than one, the 
resonance condition can be fullfilled only in a highly neutronized medium, which can be 
found in supernovae explosions. In Ref. [|1^ it was shown that the neutrinosphere appears 
for the electron fraction Fg ~ 6 x 10~^. The experimental upper bounds (^ show that 

the resonance condition can be fullfilled for Fg ~ 0.015, quite close to the Fg value at the 
neutrinosphere. Assuming there is no nonforward scattering of neutrinos |^ , the authors of 
Ref. 10 found the probability for i/g •«->■ Vr and Ve ■^^ i^r conversions in a simple Landau-Zener 
approximation pO|,|iT|], 



P = P[ye 

^ 1 — exp I — 



V^) = l-P{Ve 



1 
2 



1 

exp 



IT 



2 Ll 



= 1~P, 



LZ, 



6r 



2 Ll 



cos 20 cos 20 n 



(9) 



where Plz is the Landau-Zener crossing probability, LJ^* is the neutrino oscillation length 
in matter at resonance, Om ~ 7r/2 is the mixing angle in matter at production point (neutri- 
nosphere) and 5r = 2sin26'|(ilnFe/(ir|~g^^. The approximate equality in (|) is a conseqence of 
the small mixing angle [0) approximation. Using that result, the expression for the detected 
terrestrial flux 



42 



;i-p) 



p 



(10) 



(0p and 0p are z/g and Ur fluxes in the absence of the neutrino conversion, respectively), the 
model independent result for the probability for u^ <->■ Uj. conversion, P < 0.35 |^, and the 
density profiles for Fg from the Wilson supernova model, Valle and his collaborators found 
the result given in Eq. (|^. 



Following the procedure of Ref. [Tj], a similar analysis can be done for the three massless 
neutrinos. To analyze the terrestrial flux data, one should know only the survival probability 
of the electron antineutrino P(z/e — >■ i^e) ^^^^- The Eq. ( ]T0| ) is still valid, but 0°^ represents 



the sum of u^ and Ur fluxes. In the three-neutrino case there are two resonances: Ue ^^ 
Ufj_ resonance and z/g ^->- j?r resonance. According to the limits (^ and the Ye value at 
the neutrinosphere, the Ue ^^ v^ resonance is within the neutrinosphere. Therefore, the 
effects of this resonance do not contribute to P{i'e —>■ ^e)- Taking that into account (or 
equivalently taking the neutrinosphere as a source of neutrinos) and using the approximative 
Kuo-Pantaleone treatment for three neutrino oscillations |^T| adjusted for physical situation 
studied here, one obtains the following expression for P = 1 — P(z/e -^ u^), 

P = 1 - {\Uel\^PLZ + (1 - PLz)\Ue2\^)){\Uel\^ + \Ue2\^) ' Pesl^ 

= 1 — Plz cos'' COS 2u; — cos'' sin^ u — sin^ u (11) 

(neutrino states 1, 2 and 3 are mainly u^, Vt and z/^ flavor states, respectively; the angles 
oj and perform rotations between 1 and 2 states, and 2 and 3 states, respectively). The 
Landau-Zener crossing probability Plz can be obtained from the Plz for two-neutrino oscil- 
lations, replacing sin 26' with 2[/ei^e2 = cos^0sin2cc; in the two-neutrino Plz- In the small 



angle approximation, assumed in Ref. |14], the probability P tends to zero only if Plz is 



almost equal one. Using the result of Ref. W^ mentioned above, P < 0.35, the small angle 



approximation, and the analysis of Ref. |14], one finds limits on mixing angles u and 



sin^ 2cu < 1 X 10-^ 0^ < 0.27. (12) 

The first limit corresponds to the limit obtained in the two neutrino case. The second one is 
too weak to give limits on the Busi matrix elements. Therefore, astrophysical measurements 
give no new limits on the heavy neutrino part of the B matrix. 

The second limit has to be compared with the v^-v^ mixing angle obtained from the 
favorite interpretation of recent Super-Kamiokande results 0, Qvf,vr ~ ^/4- Obviously, 
these two results are in a slight contradiction. 



C. Theoretical limits 

If one wants to work in the perturbative regime of the theory, an additional constraint 
on the Bijsj mixings comes from the theoretical argument that the partial wave unitarity 
(perturbative unitarity) has to be satisfied. From the perturbative unitarity follows that the 
decay width of any heavy neutrino has to be smaller than a half of its mass. Written in 
terms of heavy neutrino masses and Bin-s, this condition reads |]10| 



na 



■m 



N, 



Y.\Bi^Nf<—M^ ^ ml. (13) 

j=l OiW 

rriD represents the upper value the Dirac mass may attain in the neutrino mass matrix. The 
perturbative unitarity bound (PUB) inequalities (|TBp give upper limit on a combination of 
a heavy neutrino mass rriN- and the matrix elements P^at.. Using Eq. (0), these relations 
may be combined into the limit for the lightest heavy neutrino mass 

8 



riR 



<<«j'(i+Epn' (14) 

i=2 

where {iti%J^ = 4M^/(avKSj'fi(s2')^) and pi = m]^Jmis[^. Concerning the calculation of 
BRs, the bound is very effective if the heavy neutrino masses are equal. If the heavy neutrino 
masses differ considerably, the bound is not very restrictive. Namely, if one of the heavy 
neutrino masses is smaller than m%_^/nj( , the others may acquire infinite values not followed 
by infinitely small values of the corresponding Bin. mixings. That leads to divergent BRs. 
Therefore, one has to use the original inequality (|T3|) to restrict model parameters. One 
cannot obtain closed expressions since the model has too many free parameters, but one can 
write two very rough bounds [|T0| 



\BlN^<s'i:, 

I 2Mw _ (0) ...x 

originating from Eqs. (^ and (|l^), respectively, which have to be satisfied simultaneously. 



If the heavy neutrino masses differ considerably, the bounds ([151) are better for finding upper 
bounds of BRs than the bound (0) . 



The "realistic" Bin^-s which automatically satisfy the PUB-s ([TsD and fullfill the relation 



J2iBiNiBij^, < {s'2) Kiay be obtained by putting 

BiN. = {{s'2r' + iBi%)-'rW^'- (16) 

This choice of Bi^.-s is used below to give an estimate how large error can be done in the 
evaluation of BRs using the rough upper bound procedure presented above. The -B/at. defined 

in Eq. (|16D begins to differ considerably from the value s^' for niNi ~ 100Mvy(0.1/s^'). 
Therefore, for rriNi values smaller than 2000 GeV, the -B^at.-s are determined by experimental 
upper bounds (H) and not by the theoretical PUB hmits BijJ, . 

D. Upper bound procedure for LFV amplitudes 

The equations (|I3|) and (|TBp are the basis for evaluation of the LFV amplitudes. The 

evaluation based on Eq. (|l^) gives the upper bounds on absolute values of the amplitudes 

To| , which have to be satisfied by any model with additional heavy neutrinos. It uses the 



Schwartz's inequality for the product of two vectors. It always gives larger estimates for 
an amplitude than the approach based on Eq. (|16|) . In both approaches the phases of the 
-B/AT.-s are neglected, but in a different manner. In the first approach the upper bound value 
of the amplitude is formed, while in the second the -B^at.-s are taken to be real and positive. 
Both approaches explicitely show that the very heavy neutrinos are decoupled. That is, 
they have no influence on the amplitudes of low energy LFV processes, in accord with the 
Appelquist-Carazzone theorem and its generalization P^^. 

Here, the improved version of the upper bound procedure introduced in Ref. |jTO| is given. 
The low energy LFV amplitudes may be written in terms of 



nji 



^B*i^,Bi,Nj{,Ni. 



1=1 



no no 

EK,^.K;-dJ(w„---) and EKd,K,,Mv), (17) 



where f{Ni, ■ ■ ■), f{uj, ■ ■ ■) and f{dj, ■ ■ ■) are expressions comprising the loop functions. The 
dots represent the indices not written exphcitly. Namely, the amplitudes often contain more 
than one sum over neutrino or quark flavors. Using the inequalities that can be derived from 
Schwartz's inequality, 

\^aibiCi\ <^\ai\\bi\\ci\, (18) 

i i 

\JZaAc,\ < |a||b|(c) + |a||b|(f: |q - {c)\^Y'\ (19) 

1=1 i=l 

^), one can write the following upper limits for the 



((c) = Yh=i Cijn) and definition of s 
expressions (p!7|). 



riR 



i=l 
riR 

i=i 

^R 

|EK.,v4/(rf„ 



< 



s^s'i (!(/(• ■ ■))n\ + [E(/(^- ■ ■ ■) - (/(■ • ■))Nn" 



(20) 



i=l 



Ur 



nR 

<ElKdJ|KdJ|/(c?„---)l, 



(21) 



where ( )Ar represents the average over heavy neutrinos. The inequality (|1^) gives the best 
estimate for the upper limit if the components Cj differ considerably. The inequality (|19D 
gives the better estimate of the upper bound if the components Cj are approximately equal. 
As the amplitudes /(mj, ■ ■ •) and f{dj, ■ ■ ■) depend strongly on quark masses, Eqs. (|2l|) give 
good estimates for the upper bounds. Eq. ( pOD is effective if the heavy neutrino masses are 
nearly degenerate, because most of the /(A^i, ■ ■ ■) functions depend strongly on the heavy 
neutrino masses. If one or more heavy neutrino masses differ considerably from the others. 



then Eq. (2C) may lead even to a divergent result as the heavy neutrino mass(es) tend to 
infinity. To avoid such undesirable behavior, one has to use the combination of the upper 
bounds (|18]) and (p!9|) for each set of heavy neutrino mass values in the following manner. 
First, the heavy neutrino masses are arranged in increasing order. The arranged masses are 
devided into two sets, one containing smaller masses and the other larger masses. There are 
J + 1 such partitions, where J is a number of different heavy neutrino masses. Then J + 1 
different upper bounds of the expression YJi=i B^n BvNif{.Ni, ■ ■ ■) are formed combining the 
upper bounds (|l^) and ([T9|) , 



riR 

E 

i=l 



"^Bi^^Bi'NjiNi, 



<^M'(K/(---)).l + E(/(^^. 

is 

+ T.BLBf,^Jf{N,, 



{fi---))s)r') 



'ibi 



(22) 



«6 



where X^i^ sums over the lighter heavy neutrino masses, and J^i^ over the heavier ones. 
Finally, the numerical values of the J + 1 upper bounds ( p2|) are compared and the smallest 
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of them is taken to be the upper bound value. For amphtudes containing sums over two 
(heavy neutrino and/or quark) indices, the procedure is essentially the same. Again one 
looks for the minimal upper bound value between upper bounds obtained for all possible 
partititions of heavy neutrino masses. This procedure gives convergent results for absolute 
values of the amplitudes, and it leads to the decoupling of the very heavy neutrinos. 

It should be noted that the above upper bound procedure gives upper bounds for BRs 
for neutrinoless LFV processes. Recently, lower bound limits for r lepton decays were 
found using the Super-Kamiokande data on atmospheric deficit of u^^, and interpreting it 
in terms of the best fit to these data |M . The mild GIM mechanism suppression, coming 



from logarithmic dependence on light neutrino masses, appearing in r — >■ fil~^l~ /fip^ decays, 
leads to the lower bounds of the BRs as large as ~ 10"^''. As the experimental upper limits 
on these processes are of the order of ~ 10~^ this lower limit is very welcome, because it 
strongly restricts the window for the heavy neutrino LFV effects. However, these results 
have to be taken with caution, as the standard interpretation of the Super-Kamiokande data 



is not the only one ^^, although recent papers pH] , ^ showed that the energy dependence 
of the oscillation wavelength strongly supports the standard interpretation. It should be 
noted that the used V model can easily be modified to include masses for massless neutrinos 
P(]| . The results for the neutrinoless LFV decays almost do not change if light neutrino 
masses, consistent with Super-Kamiokande measurements, are introduced. 

IV. NEW RESULTS ON LOW-ENERGY NEUTRINOLESS LFV PROCESSES 

As mentioned in Introduction, heavy meson neutrinoless LFV decays and M ^-y M 
conversion have not been studied in the V model. They are examined below. Some previous 
results for neutrinoless LFV decays are extended and/or corrected. 

A. Neutrinoless LFV decays of heavy mesons 

The LFV decays of heavy mesons were discussed in a few papers in the context of the 
leptoquark models |T^ , a flipped left-right symmetric model [|T7| and SM with an additional 
Higgs doublet [|15|. 

In these decays both lepton and quark flavor are changed. In the V model they can 
proceed only through box diagrams in which two W bosons are exchanged. The effective 
Lagrangian on the quark-lepton level reads 

^e// = T^ E E E Fl'^t" H{i - i.)i' 

Qal^il - 75)Q [5qc5,^u - SQbiS.^d + Ks)]- (23) 

/ and /' are the lepton fields, qa and Q are the light and heavy quark fields, respectively, aw 
is the weak fine-structure constant, Mw is the W boson mass and Fj^^^ is the composite 
loop function, 

riR no 
^box" = E E B*ii N^BiN.V*^Vcd^ [Hboxi^NiAdj) 
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—Hbox{^Ni,0) — Hbox{0, AdJ + HboxiO, 0)], 

nR riQ 

Kox"' = X! X! ^I'N, Bin, Vu^da Vu,b [F'boxi^N, , Kj ) 



-Fboxi^N,. 0) - Fbox{Q, A„J + Fbox{Q, 0)]. (24) 

Fbox and iffeoa; are loop functions defined in [0]. These loop functions have approximatively 
logarithmic dependence on the heavy neutrino masses. 

The dominant processes are those which have maximal value of the CKM matrix ele- 
ments, maximal LCKM matrix elements and t-quark mass in the loop function. The main 
neutrinoless LFV candidates, between the two-prong and three-prong processes studied here, 
are B^ -^ r^e^ , B~ — > K~T^e^ , B^ -^ K'^r^e^ and B^ — >• (pr^e^ . There are no interesting 
D meson candidates for two reasons. One is of dynamical origin - the quark masses involved 
in loop functions are smaller than in B-meson decays, so loop functions are much smaller. 
The only larger loop contribution coming from the b-quark is suppressed by small CKM 
matrix elements. The other is kinematical - the difference of r lepton and D meson masses 
is small. The small quark masses in loops and large t-quark width makes LFV decays of 
t-quark uninteresting from the experimental point of view. 

The matrix element of the neutrinoless LFV decay of a heavy meson H, H —>■ XW, 
contains hadronic matrix element {X\qa{0)'~f^{l — '~f5)Q{0)\H) . The corresponding matrix 
elements are usually parametrized in the following way [pl|,0] 



(0|g,(0)7^(l - 75)Q(0)|i/a(p)) = -if HP, 

{Pip')\qai0h,{l-l5)QmHa{p)) = 



ip + pX-^^^^^^<l.)F^{q^ 



+ !±_!!±q^F,{q^^ 



Np, 



{V{p',e)\qMr{^-l5)QmHa{p)) 



^-^^^^e^-^elp^p', - ^e* • g^g,Ao(g^) 



niH + my q 

{{p + p).-'^^^^%)A.{q 



niH + my \ q^ 



+i{mH + my) ( £* —q^, ) Ai{q 



N'y. (25) 



Ha is a heavy pseudoscalar meson containing light quark g^, P and V are a light pseudoscalar 
meson and a light vector meson, respectively, p and p' are 4-momenta of the heavy and light 
meson, respectively, q = p — p' is the momentum transfer, e is the polarizaton vector of the 
light vector meson, fn is the decay constant of the heavy pseudoscalar meson, Fi, F2, V , 
Aq, Ai and A2 are form factors and Np' (Ny) is a factor in front of the term containing g^ 
in the quark wave function of the P {V) meson. The q"^ dependence of the form factors is a 
consequence of long-distance (resonance) effects following from strong interactions. 

To evaluate the hadronic matrix elements of quark currents and to include the long 
distance effects, one has to express the quark currents in terms of the meson states and 
to introduce a strong-interaction Lagrangian on the meson level. Similar hadronic matrix 
elements have been extensively studied in radiative, semileptonic and nonleptonic decays of 
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heavy mesons. The combination of heavy quark effective theory (HQET) and chiral pertu- 
bation theory (CHPT) has been apphed to these decays ^3\. Here, the modification of this 
formahsm [^ , p^ -|56| is used. The authors of these papers replaced the HQET propagators 
by the full heavy-quark propagators, and introduced SU{3) symmetry breaking through 
physical masses and decay constants of light mesons. The matrix elements in that approach 
read 



(0|g„(0)7^(l - 75)Q(0)|H(p)) = -tfnP, 
{Pip')\qaiOh,il-l,)Qmnip)) = N: 



ip 



X9[Pu- 



J 7i \ c\J' 

P' ■ (i% 



Jn , nJn'* f Ni/2 

- -rPt^ + 2-—{mnmn") ' 
Jp Jp 



m 



n' 



mw 



q^ — m|^,: 



{v{p',e)\uoh'{i-i,)Qmnip)) = N^'^ 



2'/' Xgy ("'''' 



V nin 



1/2 



/. 



H'- 



m-u' 



q" — m^. 



2 e,,uape*^py'^ 



:2^/V^v 



rriH' 
m-H 



1/2 



/' 



W 



q ■ £ ^A* 



q' 



m 



1 2^/^aigvmlih*'' 



. .^1/2 1/2 PmP • ^ 



m 



n 



(26) 



7i' and 7i'* represent heavy pseudoscalar meson and heavy vector meson resonances, 
respectively, fn', fn'*, "^w, m?^/* are the corresponding decay constants and masses, 
gv (~ 6.0(2/a)^/^ with a = 2 in the case of exact vector meson dominance) is the vec- 



tor meson self- interaction coupling constant [0, g and /3 are the coupling constants in the 
even part of the strong interaction Lagrangian |^,^|5i|,^|5S|,^ , A is a coupling constant 



in the odd part of the strong interaction Lagrangian |[5T|j5^j5^ - |56| , p8| , and ai and 0^2 are 
coupling constants in the definition of weak current [|^,^. The constants g, (3 A, ai and 
a2 are free parameters which have to be determined from experimental data. 

The matrix elements of 7i — *> XIV follow from (0), (|2^) and (^61). From these matrix 
elements follow the corresponding decay rates: 



Binl -. /-/'+) 



«VK fno''^no ^^^'^{^noif^lif^D m^oimf, + mf) - (mf, - m 
2^HTnoM^ 



:pl'lqaQ\2 
box I ' 



n° 



no 



B{Ha 
B{Ha 



vi-i'+) 



c^UNpr /(W 



+m;,)2 



dt [apZpi + apbpZp2 + bpZp^] 



13-7r3 



21% 



a 



w 



{N'v 



2i%3 



\'pl'lqaQ\2 
I 60a:: I 



1 [■{m-H-mv)'^ 



ml^T-nM^ J(mi+mi,Y 



dt 



■^l'lqaQ\2 
box I ; 



O-V^Vl + byZy2 + CyZy^ 



-\-dyZy4^ + ayCyZy^ + OyCyZyQ + bydyZyj + CydyZ 



V8 



(27) 



The form factors ap, bp, ay, by cy and dy, and phase functions Zpi 
i = 1,- ■ ■ ,8 are defined in Appendix. 



1, 2, 3 and Z 



Vi 



B. Muonium antimuonium conversion 

The CC vertices in the V model have V — A structure. The effective Lagrangian for the 
M ^->- M conversion comes from the lepton box amplitude. Therefore, the structure of the 
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effective Hamiltonian density ioi M ^^ M has the same {V — A) x (V — A) form as in the 



Feinberg's and Weinberg's papers [^ 

^ = Gmm V^m7a(1 - 75)V^e^M7^(l - 75)^^6 



(28) 



in which they had elaborated the original idea of Pontecorvo |]6T[]. The constant G^m 
contains information on physics beyond SM. In the frame of the V model it comprises the 
parameters of the box amplitude for the process /i"'"e~ — > fi~e~^, which is forbidden in SM, 



G 



a 



MM 



W 



_ p/^ee/x 



WMy""" 



(29) 



box 



is a composite loop function having the following structure fTL 



riR 



(30) 



Using the expression 



for large degenerate heavy neutrino masses, one obtains the limit 
Gmm < 3.9 x 10"^,, Gp, (31) 



p^J.ee^J. 
box 



/(0.5A 



nKSl 



'm\2/ 



From the definition 



follows that the a;?,„^„ may assume only values 



"fieefj. 



where Gp is the Fermi constant and a;°gg^ - 
of the composite loop function and the limit 

smaller than 4.7 x 10~^. Keeping that in mind, the result (|3T|) has to be compared with the 
recent experimental upper bound [^,^ which improved the previous experimental result 
50, Gmm < 3.0 x lO^'^Gp- The upper bound ( |3TD is larger than the result 
estimated within SM with massive Dirac neutrinos, by comparing the 



|63| by the factor 



found by Swartz |64 



effective Hamiltonians for M-M conversion and for B^-B^ transition. Having in mind that 
the upper limit (^ was much weaker than when Swartz wrote his paper, the result obtained 
here is in fact larger than the numerical results show. The G^m w^ls also evaluated in many 
other models |6^. Depending upon the variant of the model, the value of G^m ranges from 
lO-^G^ toO.lG^. 

The conversion probability P{M -^ M) is the quantity that is measured in experiments. 
It is related to the constant G^m ^^ the following way 

5^ 
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P{M -^ M) 



2^2' 



(32) 



mm 



(33) 



where 

is a transition matrix element between the muonium and antimuonium states (a is the radius 
of muonium atom) and F^ is the total decay width of muon. 

From the point of view of SM extended by heavy neutrinos, M-M conversion is not a 
good place to search for LFV. Roughly speaking, the M-M amplitude is proportional to 
the square of the nondiagonal /x-e mixing X)j B^NiB*j^^, which is strongly constrained by the 
measurements of processes // -^ 67, /x —>■ eee and /i — >• e conversion. Amplitudes of the three 
processes depend approximatively linearly on the nondiagonal fi-e mixing. Therefore, if any 
of the experimental results of the three processes is improved by a factor a, the experimental 
result for P{M -^ M) has to be improved by the factor a^ to be competitive in finding LFV. 
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C. Extension and correction of some previous results 



In this subsection some previous results on neutrinoless LFV processes evaluated within 
the frame of the V model are extended and/or corrected. 

The decays of r lepton into three leptons were evaluated within the V model in Ref. 
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without including terms with four Bin-s. These terms were shown to dominate for large 
heavy neutrino masses in SM extended by two additional heavy Majorana neutrinos |11]. In 
that model the Bu^-s are completely determined by s^' parameters and ratio of the heavy 
neutrino masses. Here, the upper bounds of complete amplitudes are evaluated within the 
V model, and used to find the upper bounds of the corresponding BRs. 
Neutrinoless LFV decays of the Z boson were studied in Ref. 



11 in SM extended with 



heavy Majorana neutrinos. The expressions for loop functions are given in Appendix A of 
that reference, and they are correct except for terms containing 



(34) 



w _i 

tan 



w 



Aj + Xj — \z 



which should be replaced with the expression 



e{w) 



A; 



+d{-w) 



tan 






\ + A, — A; 



+ Tie{\z -\i- Xj) 



In 



\i - A,- + ^A^ 



Az - Aj - Xj - y/^W 



- in9(Xz - Xi - Xj) 



(35) 



The notation is the same as in Ref. [jrT| . The theta function in the first square bracket was 
not taken into account in the analysis in Ref. |ll|]. As it contributes only for the heavy 



neutrino masses smaller than the Z-boson mass, the numerical results given there should 
not change. For heavy neutrinos lighter than Z-boson mass, the theta function assures the 
continuity of the loop functions in heavy neutrino masses. Here, LFV decays of the Z boson 
are studied in the V model. The terms containing the matrix elements Cj^jy., that exist 
only for heavy Majorana neutrinos, are neglected. In the V model only the upper bounds 
of the Z — > //' amplitudes can be found. They are found using the formalism of the section 



IHC 



The only three neutrinoless LFV processes that give additional constraints on -Bjat-s, 
fi —>■ 67, fi —>■ eee and fi — e conversion, were examined in Ref. [0. Their analysis has 
included the " nondecoupling" effects of heavy neutrinos, has indicated that a generalization 



of Appelquist-Carazzone theorem ||4^ , |45[| is valid for the V model and has determined the 
limits on specific combinations of Bi^-s. The "proof of the generalization of the Appelquist- 
Carazzone theorem is based on an introduction of a somewhat arbitrary maximal SU(2)l- 
doublet mass term. The amplitude they present for /i — ;> e conversion does not include the 
photon exchange and box contributions, and the amplitude for /x — >• eee does not include 
the photon exchange term. These terms are included here. Moreover, in their expression 
for /i — > 3e BR, obtained in the limit of large heavy neutrino masses, one has to make 
replacements J-'e^ 



2Ten and el 



-l/2 + s'i^ 



-El, (the notation of Ref. [jT2| is used). 



The neutrinoless LFV decay of tt^ was studied in Ref. [^ in extensions of SM with 
additional Majorana and additional Dirac neutrinos. The expressions for the extension with 
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Majorana neutrinos is correct, but the expressions for the extension by Dirac neutrinos is 
not, because the terms existing only for Majorana neutrinos were kept in the aniphtude. 
The correct aniphtude is obtained neglecting the terms containing the loop function Hz- 
When this correction is made, the numerical results for the vr — >^ /xe decay become ~ 25 
times smaller. 



V. ON LOW ENERGY NEUTRIONOLESS LEV AMPLITUDES AND DECAY 

RATES 

A. Loop functions included in LEV processes 

In the lowest order of perturbation theory, amplitudes of neutrinoless LFV decays are 
built up from several building blocks (composite loop functions and tree-level functions) 
which may be denoted by the exchanged bosons, or by the type of the Feynman diagram: 
7, Z, box (box containing only leptons, leptons and u quarks, leptons and d quarks), H 
and W^W~ . All functions except the last one are combinations of loop functions and Bin-s 
n| , pil|J5^ . W^W~ function is a tree-level function and it is strongly suppressed compared 



to the others ||6^. 7, Z, box and H functions comprise two-fermion currents. In the 7, 
Z and H functions only one of the fermion currents changes flavor, while in box functions 
flavors may be changed in both fermion currents. The classification of the neutrinoless LFV 
decays, given in Table I, is made according to the Feynman diagrams they contain and 
the approximations (physics) one has to use in finding the corresponding amplitudes. The 
references cited in Table I refer only to the calculations of LFV processes in the extensions 
of SM by additional heavy neutrinos. 

If the heavy neutrino masses are larger than a few hundred GeV , the expressions for neu- 
trinoless LFV decays simplify considerably. All amplitudes can approximately be expressed 
in terms of four combinations of masses and Bin^-s, 

Aw = "^Bij^Mi'N,, 

Bui = Y^ Bij^Mi'N, In A at,, 

r - sr n* r* R ^N,>^N, , Atv, 

Wr - 2^ ^lN,^N,N,^l'Nr^ T ^^\ ' 

T^Whh = o H ^in.BInABi'n.Bi^n, + Bi^n,Bun,)t — / ^'^T^' (^^) 

^ 7V,7V, -^N, - Aat^. An, 

where Xn, = i^'nJi^w- The building blocks mentioned above, expressed in terms of combi- 
nations iWi). read 



G^ ~ -Aw, 

Fz ~ ~2^"' ~ 2^"'' 
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rpll'hh r^ 
box ^^ 



-pU'UaUi, 

box 






+ 



A. 



KT.feKtfe 



i3/r 



box 



- ^dadt 



Ui^C^t 



?,i- 



4(1 -A,) 



3 A„^ -A^^ + 8A^^ - 4A„^ 



Aw 



X. 



+ 



4(1 - A„, 



lnA.jK.d.K:. 



Au> 



+ 



A,, 



E-— IT/ T/* 
ni=c,t 



i3«. 



7:11' 

H 



G 



ir 

H 

no 



A 



H 



7:Au' + -^Eiv + -A 



Fw+w-- {T.^u,d^'ld)jAw 



(37) 



where A^; = m^^/m^, x = b,t, H. 

For the important case of degenerate (Aat. = Aat) and large heavy neutrino masses the 
functions (|36D can be written in terms of parameters s'2 and x°;/, 



Aw = sis I Xii,, 



N, 



riG 



c.' = 3^2'E(^L)'44.Aiv, 



i=l 



'^ll'hh = ^^L^Ln'-H^i^ll'^hh + ^lh^Vh)^N- 



(38) 



It is convenient to introduce four combinations of Bin-s, heavy neutrino masses, A^^"^ and 
upper bound values for s^ parameters (^), denoted by s^: 

xw = Aw{s2s'2') , 

zw = Bw{srs1'\nX^^^^y\ 

?/.'=C,.(5M"E(Sl')^A^^^)", 



yii'hh 



i=l 






(39) 



Any of these combinations is always smaller than one. 

Here, few comments are in order. First, it is obvious that I'Dwi^iil < \Cw\ (the relation is 
also valid for large, nondegenerate heavy neutrino masses). Second, for degenerate neutrino 
masses, the function Cw becomes larger than the functions Aw and Bw if 



Atv ~ 



JIMs'l)' 



and Aat ~ 



In A 



N 



n^iis'i:r 



(40) 



respectively. The dominance of the functions with quadratic mass dependence of the ampli- 
tude leads to the transient, so called "nondecoupling behaviour" of amplitudes. As explained 
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in section [III D| , decoupling follows from PUB inequalities ([T3[). A typical mass value for 
which the quadratic terms become larger than the logarithmic terms is rriN ~ 1500 GeV 
for s^ values of the order of the present experimantal bounds (^). Third, at the maximal 
\n value permitted by the PUB (A^^"^), the function Cw depends essentially only on two 
diagonal mixing parameters, s'2 and s/ , 






(^wlli[s£Y aw aw 



Therefore, at ttiat = m^^ all amplitudes depend essentially only on s'2 and s^'- If both the 
logarithmic and quadratic mass terms are present in LFV amplitude, at m^^ logarithmic 
terms contribute up to ~ 10% of the total amplitude. Fourth, if the t quark contribution 
is multiplied by small CKM matrix elements, box amplitudes may have large contribution 
coming from c quark in the loop expressions. For instance, in the processes r -^ eP^/^P^ 
c quark contribution to the amplitude is ~ 13%. Fifth, the processes containing only the 
function Aw are most suitable for obtaining new information on Bi^i parameters, because 
they are almost independent of heavy neutrino masses. Sixth, for degenerate heavy neutrinos 
the dependence of LFV amplitudes on LCKM matrix elements appears only through six sums 
J2iBiNiBl,j^., I 7^ I' and J2i IBin.I"^ (diagonal and nondiagonal mixings). Writing the sums 
in terms of s^-s and x^,-s, one can easily show that if some LFV amplitude tends to zero 
for s^ — > 0, then the amplitude tends to zero for x^, ^ 0, I j^ I', too. (Stricktly speaking 
reduction of a parameter s'2 by factor a is equivalent to the reduction x'li — > a'^x^ and 
x'li, -^ ctx'ii,, I ^ I', but, by definition, a;°; = 1.) This analysis shows that LFV amplitudes 
may be reduced without changing the diagonal mixing parameters s^' . It also indicates that 
the absolute values of LFV amplitudes may attain any value between zero and the upper 
bound value. Therefore, LFV processes are unsuitable for finding the limits on the diagonal 
mixing parameters s'2- 

B. Approximative expressions for BRs in the large- mass limit and relations between 

them 

Keeping only the leading terms in the large-mass limit of heavy neutrinos, the expressions 
for BRs of neutrinoless LFV decays may be expressed in terms of the functions (pBI) . In the 
following these expressions are listed. The definitions of unknown quantities are given below 
the list. 

3 2 5 

R(l -^ I'^A ~ '^W^W ^l I A |2. /.9^ 

B{l^ll) ^ -^I^J^J^i^'l ^ (42) 



4 5 1 

Oi-w ^l (\^ 0/1 o„2 \n |2 , ^\A„2 n |2 , 



Bil- ^ /'-/-/+,/' = h = k) ^ 3x2i6^3M^rA '^"''' " ^^^ " ^'w)Cii'\' + ^l^sl^Cw 



aw mr 



B(r ^ ri^itk # i:h) « ^^J^^ lA,.,.,/; (43) 
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M, 
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3 X 2^clr Tz 



r |2. 



(44) 
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m. 
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\c P- 



(45) 
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(51) 
(52) 
(53) 



P(r^/PiP2,cqf) 
P(r^/PiP2,ncqf) 
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!{Z+^2? ^^ " I ^V-o Pw(9)"f oCyop^p, I 



216^^3 



box,sd\2 
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box 
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BiHa ^Vl I } IJ^box I TTTZZT-FT- . ^ ., ^t [ayZyi + byZv2 

+CyZv3 + dyZvi + ayCyZv^ + hvCyZv^ + hydyZyi + CyrfyZyg]; (59) 



air ,_„.rf„,o 1 /•(ms-ms')^ 



+A,{fm) + A,{gigs) + A,{gl)]. (60) 



All expressions are written in terms of products of diniensionless factors. For the expressions 
containing the dominant term Cui, the error one makes by keeping only the dominant term 

is of the order ~ 20%, because the term Cui is always accompanied with the Bw term 
giving ~ 10% contribution to the amplitude at m^^ . Following abbreviations are used: 
sw = sm9w, cw = cos^vi/ [6^ is Weinberg's angle), sp = sin 9p, cp = cos9p [Op is the 
mixing angle for psudoscalar nonet states) and sy = sin 6'p, cy = cos 6'p {Oy is the mixing 
angle for vector nonet states). In Eq. ( |45| ) aem = 1/137 is the fine structure constant, Z 
is atomic number (for |^Tz Z = 22, N = A - Z = 2Q), Z^jf = 17.6 0-|7| is the effective 



atomic number of Ti [^], F{—m^^ = 0.54 is its nuclear form factor 0,^] at momentum 
transfer g^ ^ — "^^ ^^, Qw = Z{l — 4:s'yy) — N is the coherent nuclear charge associated with 
coupling of Z boson to nucleus [^ and T capture is the capture rate for negative muons on Ti 
||71| , [76|J77| . In Eqs. (|47|-|53D fp and fn are decay constants of light and heavy pseudoscalar 
mesons respectively, and 7yo are constants defining the decay constants for light vector 
mesons, fv = my j^y. The normalizations used here are 

A^p{x)=%fpd^P{x), 

.,,/t/ N I ^^y(a;) for light vector mesons, ir.-,\ 

Vy\X) = \ 'IV ^ > ^ (61) 

I fymyV'^{x) for heavy vector mesons, 

where Ap and Vy are the axial vector (vector current) with the same quark content as 
corresponding pseudoscalar meson P(x) and vector meson V{x) fields, respectively, apo, 
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Oiyo, ctpo' ° '' and ayo' ° '' are constants defining the meson content in vector and axial- 
vector quark currents contained in quark combinations in the Z- and box- amphtude. They 
are defined in Table II. The mass of the lighter lepton is neglected in the expressions (^^ 
1531) . As the composite loop function F^^^"'^^ contains two terms of approximatively equal 
magnitude, for brevity the expresions of these BRs are not written in terms of the functions 
(P5|) . The expressions (^3|-^) contain three body phase integrals. The phase functions may 



be found in the following references: a, /?, 7 and l in Ref. [£6|; A_(.+ , y4_|__ and A in Ref. 

n| ; Ai, A2, A3, A4 and A^ in Ref. [0. In Eq. (^^, Cyop^p^ = S'pOir+Tr-Cvopjp^ are constants 
stemming from the gauged chiral U{?>)l x f/(3)ij/f/(3)v' Lagrangian |^ (e.g. CpO^+^- = 1), 
and 

A(g) = ^^^iZilVom^ ^ —pVr-{^) (62) 

is the Breit-Wigner propagator for a vector meson V^ multiplied by slightly modified expres- 
sion rriy/'^Y. The modification of the expression rriy/'^y is made to obtain PBW°^^i^) — 1 



,[78|-p0|. The constant gpOT,+.„- is equal to the p self couplig constant gy from section |IVA 



In Eq. ( pBD Mjjp^p^ are mass parameters contained in the effective Higgs-meson Lagrangian 



-HMM 



_9W_ 

4M, 



w 



m^ ((7rO)2 + 27T+7T-) + 2m%^K+K- + 2m^oK°i?° + ^^++<° + ^' ^2 



2mi+ + 2mio - m^ 2 , 23/2(^2m2 - m|.+ -m 
\ i^ -^ r^g H — ^i ^^^VsVi 



2 -1 

X0> 



(63) 



obtained by comparing the quark mass Lagrangian with the corresponding term in the 
chiral Lagrangian ^ (e.g. M^^+^_ = 2ml = 2(2m2+ + mlo)/3). In Eq. (|3) 
ck*okwtt = '^c^*o_Ss7r + ^c/f*o^g^ (-K'vK = ai^s + bKs is a weak kaon eigenstate, and Ks 
is a mass eigenstate). In Eq. (|60D /i, (/i and (73 are baryon form factors. The other baryon 
form factors do not contribute, because they belong to the second class currents, or give a 
contribution proportional to the difference of baryon masses. The form factors f\ and gi 
can be defined in terms of two SU{3) Clebsh-Gordan coeficients and two reduced matrix el- 
ements corresponding to symmetric and antisymmetric octet SU{3) representations. These 
reduced matrix elements are almost independent of momentum transfer and are usually 
identified with their value at zero momentum transfer, D and F. The functions gi and g^ 



are not independent, but correlated through the Goldberger-Treiman relation [10 



All approximative expressions for the BRs (B2|-Bn), valid in the large mass limit, except 



Bi^l" -^ l'^lil2,h = h ^ V) and B{1~ -^ l'~lil2,l' = h = h), depend only on one of the 
functions (PB[). In the following, the smaller of two dominant functions T>iiq-^i^ will be ne- 
glected in the two exceptional expressions. The maximal error one makes in the evaluation 
of BRs of the exceptional processes is ~ 40%. With those approximations, the ratios of BRs 
having the same dominant function ( pG]) become independent of the l^-model parameters: 

a) BRs with F^" (C^^): 

i?(/iTi -^ eTi) : B{fi -^ ee'e^) : B{Z -^ /i^e^) : fi(7r° -^ p^e^) : B{t] -^ /i^e^) 
= 1 : 5.60 X 10~2 : 3.77 x 10"^ : 6.05 x 10~^° : 1.69 x 10"^\ (64) 
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b) BRs with Fp (Cri): 

B{Z^T^l^) : B(r^/7r°) : B{t ^ Ip^) : B{t ^ In+n-) : B{t ^ Icj)) 
: BiyT^irt) : B{t^II^I+) : B{t^IK+K-) : B{t ^ IK°K°) : B{t ^ Ir]') 
: B{T^lr]) : B{t ^ luo) : B{t ^ lr]r]) : fi(r ^ /ttV) 
= 1 : 3.40x10"^ : 3.17 x 10^^ : 2.83x10"^ : 2.81 x 10~^ 



: 2.64 X 10"^ : 1.64 x 10"^ : 1.20 x 10"^ : 7.43 x 10"^ : 6.15 x 10"^ 




: 4.72 X 10-2 : 8.78 x 10"^ : 4.34 x 10"^^ : 5.50 x 10-^^ 




c) BRs with F,^''' 




B{KL^p^e^) : 5(7^+ ^ vr+^u^e^) : B(S+^p^^e±) : 5(S° ^ A^ 


^e±) 


: B{K ^ n^^e^) : B{Z- -^ ^-^i^e^) : 5(S° -^ SV^e±) : 5(S° -^ 


njl^e^] 


= 1 : 3.01 X 10^2 : 1.30 x 10"^ : 1.21 x 10^^ 




: 8.66 X 10"^ : 6.40 x 10"^ : 4.07 x 10"^ : 6.31 x 10-^^ 




d) BRs with FC'"" 




B{B~ -^ 7i-p^e^) : 5(fi° -^ fi^e^) = 1 : 3.76 x 10-^ 




e) BRs with F/^,f 




B{B- -^ K*-fi^e^) : 5(5° -^ K*y^e^) : 5(5° -^ r]'fi^e^) : B{B^^ - 


^ (pjji^e 


: B{B- -^ K'fi^e^) : 5(5° ^ i?°/i^e=^) : 5(5° ^ ^/i^e^) : 5(5° - 


-^/i^e± 



(65) 



(66) 



(67) 



= 1 : 9.34x10"^ : 8.83x10"^ : 8.57x10"^ 

: 7.92 X 10"^ : 7.47 x 10"^ : 3.31 x 10"^ : 4.93 x 10-^ (68) 

f ) BRs with Ff^jf 

B{t -^ en^K') : B{t -^ eK*°) : B{t -^ eK°) = 1 : 7.32 x 10"^ : 2.99 x 10"\ (69) 



g) BRs with F;^y 



B{B- -^ TT^r^e^) : B{B^ -^ r^e^) = 1 : 1.14 x 10"\ (70) 



h) BRs with F^,^^' 



box 



B{B- -^ K*-T^e^) : B{B'' -^ K*\^e^) : B{BI -^ 0r^e^) : BiB^^, -^ ri'r^e^) 
: B{B- -^ K-T^e^) : B{B^ -^ K^r^e^) : BiB'^^ -^ r]T^e^) : BiB'^^ -^ r^e^) 
= 1 : 9.37x10"^ : 8.10x10"^ : 6.57x10"^ 
: 6.54 X 10"^ : 6.16 x 10"^ : 2.76 x 10"^ : 1.63 x 10"^ (71) 

For each group of BRs the BRs are hned up in the descending order. For instance, fx -^ e 
conversion is the most suitable for finding LFV in the group containing the composite loop 
function F^^. The position in the group depends on the coupling constants, phase factors 
and the total decay rate of the decaying particle. For instance, the BRs for r -^ IP1P2 
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processes containing Z-boson amplitude are ~ 10^^ times larger than BRs of the processes 
T -^ IP1P2 containing only Higgs amplitude, because of the small Higgs-meson couplings, 
although the dominant composite loop functions are essentially the same. The ratios ( [6^ - |7T[ ) 
are given for measured processes and LFV processes that have not been studied in models 
with additional heavy neutrinos before. The / — > /'7 decays are not included in the above 
ratios, because each process / -^ l''~f forms a group for itself, depending only on the function 
All'. The numerical results for the ratios of BRs agree quite well with the exact ratios for 
degenerate heavy neutrino masses obtained at mjv = rn^^ . That allows one to consider 
only one of the decays of each group when comparing theoretical and experimental results. 
Besides the ratios of BRs having the same dominant function (0), it is usefull to have 
relations that relate BRs of different groups of decays. These relations generally depend on 
the matrix elements of B matrix and CKM matrix. For instance, 

I TpTe\2 I T7TM12 I 77/^6 12 

- FZ I • \^Z I ■ \^Z \ ^ 

B{t -^ eP°) : B{t -^ fiP^) : 5(P° -^ e^/i^) 



\-pTeds\2 . \rpTflds\2 . 1.45xl03g-l "^pO (^""^A^Mpo) | jplJ.edsi2 fr^-r -nnrif 



B{t -^ 67) : B{t -^ fi'j) : B{^ -^ 67) 

I R* R |2 . I r)* r) \2\2 . r col R* R |2 

— l^rN-^eNl ■ \-DrN-^nN\ ) ■ 3.0-3 |-D^jv-DeAf| , 

S(r^ee^e^) : B{t ^ fifi^fj.^) : P(/i ^ ee^e^) 
^ |F|"|2 : |F7|2 : 5.63|Ff p, 
B{T^eK+TT~) : B{t ^ nK+n-) : B{K- ^ n- fi^e^) 

|7-ireds|2 . n noQl T7TlJ-(is\2 . c o o I l?Meds 1 2 

- Ffeox I • ^-^^^l-^box I • 0-^^Ffeoz I 5 

5(5. -/?/f) : 5(5, -/^/f) = \Flt'^'\' : iF^f*? g.,g,=«,rf,.. (72) 

BRs of processes having only the logarithmic dependence on mass are several orders 
of magnitude smaller than BRs containing the quadratic mass dependent terms. In the 
processes containing quarks in the final state, the presence of small CKM matrix elements 
additionally reinforces this difference. E.g., at m^^ 



^ lo-^ 

B{t- -^ e~K^) : B{t- -^ e^vr") ^ B{t- -^ e^Tr+iT") : B{t- -^ e"7r"7r+) ~ 10~^ (73) 



B{l^l'-f) : Bil^l'hh) ~ 10 , 

n 



Between the LFV decays having the box contribution only, the 5-meson decays have the 
largest CKM matrix elements. For that reason they might be the most suitable box- 
dominated processes for finding LFV in the future 5-factories. 

C. Numerical results, comparison w^ith experiment and discussion 

In this subsection the experimental upper bounds for the measured neutrinoless LFV 
BRs are compared with the theoretical upper bounds obtained in the V model. For some 
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interesting unmeasured processes, the theoretical upper bounds are given, too. The results 
are discussed. The limit on the nondiagonal /i-e mixing is updated. The decoupling of very 
heavy neutrinos is shown explicitely. The possible error one can make using the upper bound 
procedure given in the section [1IID| is estimated. 

Theoretical results depend on the ^- model parameters: "diagonal" mixings s'2, phases of 
BiN-s and heavy neutrino masses. The parameters s^ must satisfy the experimental upper 
bounds (^, the heavy neutrinos are bound by the PUB inequalities (|TB|) and (|T^, while the 
phases of Bin matrices are undetermined. The numerical results are largest for degenerate 
neutrino masses at maximal values of s'2 parameters and maximal neutrino mass, m^^ . 
For degenerate heavy neutrino masses the phase dependence of Busf matrices is contained in 
the parameters xw. 

The numerical values for BRs and Gmm depend on a number of " SM" particle properties, 
too: decay rates of particles, masses of the particles included in the decays, CKM matrix 
elements, decay constants of mesons, quark masses included in loops, mixing angles, various 
couplings, etc. Almost all these quantities are taken from Ref. [^, or derived from the data 



given there. For instance, masses of the m, (i, s, c and h quarks are taken to be equal to the 
average of the upper and lower bound values. The CKM matrix elements are derived in the 
same way. The t quark mass is set to be equal to the experimental value obtained from the 
direct observations of t quark. For pseudoscalar meson decay constants of light mesons we 
took the values partly from Ref. |^ and partly from Ref. [^ , 



f^+ = 130.7 MeV, fK+ = 159.8 MeV, f\o = 119 MeV, 

f^ = 131 MeV and /^/ = 118 MeV. (74) 

Due to the isospin symmetry, /^o = /^o = /i^±. The constants 7^/0, defining the decay 
constants of light vector mesons, are extracted from the V^ — >■ e^e~ decay rates 

7pO = 2.518, 7,^ = 2.933, 7^ = 3.116, (75) 

or estimated using the SU(3) octet symmetry, 7/^*0 = 7p0. For all decay constants of D and 
D* mesons, the conservative value 200 MeV is taken. The decay constants of B and B* 
mesons are derived using the scaling law for decay constants derived from HQET, 

fn ~ m:^^^^- (76) 

The weak fine-structure constant is defined as aw = aem/sin^6'iy, with cos9\y = M^/Mz- 
The p — TT — vr coupling constant (which is equal to the p meson self coupling constant) 
is derived from the p — *> 27r coupling width. Other vector-meson-pseudoscalar-meson cou- 
plings of light mesons are fixed by one of the chiral model Lagrangians |l57| , |66| . The mixing 
of the vector-meson nonet states is determined from the quadratic Gell-Mann-Okubo mass 
formula, By = 39.3°. The mixing of the pseudoscalar-meson nonet states is extracted from 
the e+e" -^ e+e-77* -^ e+e-(P -^ 77) experiments [g, ^p = -23°. The only "SM" 
parameters that are not firmly established are "HQET+CHPT" parameters describing the 
semileptonic LFV decays of the 5-mesons, g, (3, A, ai and 0^2 (see section |1VA| ). The 



corresponding parameters for D-mesons have been determined by fitting the theory to the 
experimental values of the semileptonic decays of D mesons [^p4| , p5| . The 5-meson pa- 
rameters A, ai and 02 may be derived from the D-meson parameters from the scaling laws 
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for the vector and axial current |5T[]. The parameter g is independent of a heavy quark 



mass, and the value of parameter j3 is consistent with zero. The best 5-meson parameters 



obtained using the above procedure are p5 



^ = 0.2, /5 = 0, A = -QMGeV-^, 
ai = -0.13 GeV^/\ aa = -0.36 GeV^^^. (77) 

That way, all parameters are defined. 

For measured processes, the experimental and theoretical upper bounds of the exact BRs 
are compared in Table Ilia. For some interesting processes that have not been measured, 
the theoretical upper bounds are given in Table Illb. In the both tables, the numerical part 
of the theoretical results is evaluated for degenerate heavy neutrino masses and the maximal 
heavy neutrino mass permited by PUB, maximal s^ values and neglecting the Bin phases. 
The factors xw, yw and yii'ij^i2 and zw describe the deviation of BRs from these values, when 
the model parameters assume other values. The factors yw and ywiih give only the behavior 
of the dominant, m^-dependent term, on the model parameters. For m^f values for which 
the terms quadratic in Bij^f matrices begin to dominate {rriN ~ 1000 — 1500 GeV), the zw 
terms begin to dominate. 

Comparing the theoretical upper bounds for the processes of the same type with different 
leptons in the initial and final state, one can see that they are often comparable in magnitude. 
For instance, upper bounds for BRs of the processes / — >• I'j, I -^ I'hh and Z -^ W are 
of the order ~ 10^, ~ 10~^ and ~ 10~^, respectively. For that reason, the muon LFV 
processes which have been measured with the greatest precision, are the most attractive 
for finding LFV. A process with weaker experimental bounds may be interesting only if the 
parameter(s) x^i, for that process is (are) large. 

If, for a specific process, the theoretical upper bound is larger than the experimental one, 
then the process gives the better bound on a specific combination of -B/at-s than the limit 
(I). The processes for which this ratio is larger than one are /i — *> 67, /i -^ eee, /iTi -^ eTi, 
T -^ ep^, T -^ evr+TT^ and Z -^ er. For the last three processes the ratio is very close to 
one. As their amplitudes are dominated by m'j^ part of the amplitude, the new limits on 
Bin combinatons contain m% mass dependence, too, and therefore are uninteresting. For 
the first three processes the ratio is much larger than one, and they do give new limits on 



specific combinations of Bi^-s as shown in Ref. |]T2[. Since that paper was published, the 
limits on /i — > 67 and fiTi —^ eTi improved by factors 1.3 (4.1 [^]) and 7, respectively. The 
first of them gives a new limit on nondiagonal fi-e mixing, 

no 

J2 B^nMn, < 2.15 X 10-^ (1.19 X 10-^). (78) 

i=l 

To obtain the limit on the nondiagonal r-e and r-fi mixings, the present experimental 
sensitivities of r — > / decays should improve by two orders of magnitude. It is interesting 
that the pTi -^ eTi conversion also gives very good mass independent limit on the sum 
J2^=i B^N.Bl^^. Namely, pTi —>■ eTi amplitude contains mass independent part coming from 
the photon exchange. If the terms in the fiTi —>■ eTi amplitude do not cancel completely, 
one can make an estimate of the sum by attributing the whole amplitude to the photon 
excahnge part of the amplitude. That way one can only make a worse estimate of the sum. 
The limit one obtains that way reads 
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no 

£s^^,i?:jv, <3.93xl0-^ (79) 

For all processes whose amplitudes comprise only the box amplitude, the theoretical 
upper bounds are several orders of magnitude smaller than the experimental upper bounds. 
For the Kl -^ e^/i^ decay the ratio of theoretical and experimental upper bound is largest, 
1.58 X 10~^a;^g. As the present experimental limit is 2 x 10~^^ [^, its significant improvement 
cannot be expected. Although the experimental upper bounds for semileptonic LFV B- 
meson processes are weak, the corresponding theoretical upper bounds are of the order 
~ 10^^. Therefore, 5-meson decays are interesting for finding LFV decays in the near 
future. 

The recent Super-Kamiokande experiment shows there is a large mixing between i/^ and 
some other light neutrino, very probably Ur. If the additional heavy neutrinos exist, this 
might suggest a large "angle" parameter x° . Therefore, the Super-Kamiokande result might 
be a sign to search for LFV among processes with tauon and muon in the final (and initial) 
state. 

To estimate how large error one can make using the upper bound procedure from section 
1IID| , the BRs for the processes fiTi -^ eTi, Z — >• ff^r^, K^ — > e'^/i^ and B~ -^ K*~ jj^t"^ 



are evaluated using both the upper bound procedure and the "realistic" -B^at-s (|T6|). These 
processes are chosen because they have the maximal BR within the group of processes 
with the same dominant composite loop function. The first two of these processes contain 
Fj' function and the last two contain F^^^" * function only. The BRs are evaluated for 
degenerate heavy neutrino masses and two sets of s2 and x^i, parameters for which the 
maximal theoretical value for B[nTi — > eTi) is equal to the present experimental upper 
bound. The firts set is obtained from the "maximal set" (s^'-s from Eq. ^ and all xw = 0) 
by replacing the maximal value for (s^")^ with the value [s'^lY = 4.29 x 10"^'^ = 7.1 x 
10~^ X (2.459 X 10"'^)^. The second set is obtained from the "maximal set" by putting 
Xfj.e = Xre = 2.459 X 10""^. The first set is used together with upper bound procedure and 
with "realistic" Bin matrix elements. The second can be applied only within upper bound 
procedure, because the procedure with "realistic" Bin matrix elements has fixed xw values. 
In all calculations x^r is kept to be equal one in accord with the Super-Kamiokande results. 
The BRs are presented in Fig. 1. as functions of heavy neutrino mass. The figures illustrate 
the following properties of the BRs. First, for all m^v values, the upper bound procedure gives 
larger value than the "realistic" Bin-s. Second, while the BRs evaluated in the upper bound 
procedure increase in the whole region of rriN values permitted by PUB, the BRs evaluated 
with the "realistic" Bi^-s may have a maximum below the m^^ . The maximum is a 
consequence of the mass dependence of the "realistic" Bi^-s. All BRs of processes with box- 
amplitude only have the maximum, but it can appear in the BRs having the Z-amplitude, 
too. Third, the by reduction of xw-s one obtains the result which are numerically equivalent 
to the results obtained by reduction of s^ parameters. Fourth, a strong cancelation of the 
amplitude terms may appear in the BRs evaluated with the "realistic" Bin-s, as in the case 
of fiTi -^ eTi. Fifth, the error one can make in the evaluation of the maximum of BRs using 

the upper bound procedure is ~ 10 for processes with the 6oa;-amplitude only, and ~ 100 
for processes with Z-amplitudes. The flat behavior of Z — > Pl'^ at rriN ~ 100 GeV{^ uiz) 
is a consequence of treshold effects. 
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As shown in the section [III Q , all heavy neutrino masses, except one, can assume any 
value between zero and infinity. BRs should not assume values larger than one in the 
whole parameter space permitted by the model. The illustration of convergence and of good 
behavior of branching ratios evaluated using the upper bound procedure and "realistic" Bi^- 
s is given in Fig. 2. Branching ratios are evaluated keeping two masses equal, while the third 
one is assumeed to take very large variable values values. In Fig 2. the branching ratios 
for the same processes as in Fig 1. are given, but here as a function of ratio of the large 
mass {rrt]^^) and mass which is kept constant {mj^^ = rrijvs)- Graphs in Fig 2. show that the 
very heavy neutrinos decouple, and therefore, that the nondecoupling of heavy neutrinos is 
only a transient effect. For the upper bound procedure, the decoupling of the very heavy 
neutrino(s) manifests as the equality of BR values for degenerate heavy neutrinos and when 
some of masses tend to infinity, while for "realistic" Bi^-s BR-s reduce in magnitude. Figure 
2 also illustrates that the upper bound procedure is very crude in the transient region where 
the upper bounds (|T3|) and upper bound ( |T^ ) are almost equally effective as the second 
(p!5|) bound. To show that with the proper choice of the parameters experimental limits are 
always satisfied, in the first panel of Fig. 2 the additional BR curve is added, evaluated in the 
upper bound procedure for parameters for which the maximal BR value is smaller than the 
present experimental upper bound for R{fiTi —>■ eTi). Only top of the curve is seen in the 
figure. The curves obtained using the "realistic" -B/at-s are much smoother than the curves 
following from the upper bound procedure. Therefore, for nondegenerate heavy neutrinos 
good knowledge of the Bn^ matrix elements is necessary to obtain reasonable estimate of 
the BR values. 



VI. CONCLUSIONS 

All low energy neutrinoless LFV processes are studied in an extension of SM by heavy 
SU{2)l X f/(l) singlet Dirac neutrinos. The structure of amphtudes and relations between 
BRs are carefully analized. It is shown that, in principle, the neutrinoless LFV decays cannot 
give new limits on the " diagonal" mixings s'l . The approximative expressions for all BRs are 
listed, keeping only the dominant terms of the corresponding amplitudes in the large heavy- 
neutrino mass limit. The approximative BRs are compared within the groups of processes 
with the same dominant composite loop function, and within each group the experimentally 
most interesting process are found: fiTi -^ eTi, Z -^ r^l^, K^ -^ e^/i^, B^ -^ K*^fi^e^, 
T -^ eiT'^K' and B^ —>■ K*~T^e^. The upper bounds of exact BRs are evaluated using 
the improved version of the upper bound procedure found in the previous publication [[10 



The results are compared with present experimental upper bounds. For maximal values 
of model parameters, only six processes have the theoretical upper bounds larger than the 
experimental ones: fi -^ 67, /i -^ eee, fiTi -^ eTi, t -^ ep^, t -^ eTr+Tr" and Z -^ er. For 
these processes new limits on combinations of Bi]\f matrices are obtained. The first three 



have been studied before ||12[ and they give a new limit on the nondiagonal fx-e mixing. The 
limit is updated here. For the last three, the ratio of the theoretical and experimental upper 
bounds are very close to one and the obtained limit is mass dependent. Therefore, it is 
not useful. Two orders od magitude improvement of experimental sensitivities is needed to 
obtain mass independent limits on the nondiagonal r-e and r-fi mixings from r —>■ I'-f decays. 
Concerning the processes with the 6ox-amplitude only, the K^ -^ e'^fi^ decay has the best 
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ratio of theoretical to experimental upper bound. Nevertheless, neutrinoless LFV 5-meson 
decays have BRs of the order ~ 10~^, what makes them interesting for finding LFV in future 
experiments. If the structure of the massless part of the B matrix is as suggested by the 
Super-Kamiokande experiment, one may expect that in the future the processes containing 
r and fi leptons in the final (and initial) state will be most interesting for finding LFV. 
Besides BRs for the low-energy neutrinoless LFV decays, the constant characteristic for the 
muonium-antimuonium conversion G^m is evaluated. The result obtained is too small to 
be interesting experimentally. 

All above results depend only on the gauge structure of the model used and masses of 
heavy neutrinos. The results do not change if the massless neutrinos are replaced with the 
light neutrinos satisfying present experimental limits. A comment on extraction of heavy 
neutrino mixings from astrophysical observations is given. Following the V^-model assump- 
tion of massless "light" neutrinos, an analysis of oscillations of three massless neutrinos in 
the supernovae is done. The analysis gives the limits on mixings in the massless neutrino 
sector that are in a slight contradiction with the Super-Kamiokande results. 
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APPENDIX A: 

The form factors ap and bp and ay, by, cy and dy follow directly from matrix elements 
of corresponding hadronic currents ([2B|). They read 



Jp Jp m^" q - m^ 



w 



fn'* L P2-q \ (mim^,.)^ /^ 
fp \ m^,. / q^ — ni^, 



fep = 2^:^(i + ^ r'7'^^'V (Ai) 



and 



ay = 2V^A,,/; ''"r'""7'; 

q^ — mij^,* 

oy — -z ' pgvJn — o 2 ' 

q^ — m^„ 

cy = -2"^aigy{m'^f/\ 

dy = 2^/'a2 , ■ (A2) 

mf 
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The phase functions Zpj, i = 1, 2, 3 and Zyi, i = 1, ■ ■ ■ , 8 in the square bracket expressions 
in (M) read 



Zpi = / " dsi-i [2pi ■ pspi -p^-mlps- p4\ , 

^„max 
I 13 

Zp2 = / dSi3 [2{pi ■ P3P2 -Pi+Pl- PiP2 ■P3-Pl- P2P3 ■ P4)], 

gtnax 

Zp3 = (is 13 [2p2 ■ P3P2 -PA-mlp^- P4], 

for H -^ PW decays and 

dsi3 [Pi ■ PsPi ■ P2P2 ■Pi + P2- PsPi ■ P4P1 ■P2-mlp2- P3P2 ■ Pi 

min 
-13 

-mlpi -p-ipi -pi], 



(A3) 



^V2 

Zv3 
ZvA 



ds 



13 



- g^ps ■qpi-q + -q^ps ■ P4 H 2 (^2 ■ g)^ ( P3 ■ qp4 ■ q 

2 777-2 



■^q P3-Pi 



ds 



ds 



13 



13 



^ ( 1 2 

P3 ■ P4 H 2 U"2 ■ P3P2 ■ Pa - -m^P3 ■ Pi 

- q^Pi ■ P3Pi ■ Pi + 7:^Wp3 ■ P4 H 2 (^2 ■ qf { Pi ■ P3P1 ■ Pi 

2 mi * 



-l^'m^P3-Pi 



/„max 
13 
dsis [2pi ■ P3P2 ■ Pi - 2pi ■ P4P2 ■ P3], 
min 



'V6 



ds 



13 



2p3 -qpi-q + q P3-Pi-\ 0^2 ■ q{P2 ■ P3Pi ■ q 

mi 



+P2 ■ PiP3 ■q-P3- PiP2 ■ q) 



jy-j 



ds 



13 



g^Ps ■ qp\ ■ Pi - q^Pi ■ qpi ■ P3 + q^pi ■ qps ■ Pi 



H — 2 (P2 ■ qy{p3 ■ qpi -Pi + Pi- qpi ■P3-pi- qp3 ■ Pi) 

mi 



'V& 



ds 



13 



- ps ■ qpi -Pa- Pi- qpi ■ P3 + Pi ■ qP3 ■ Pi 



+ rTP2 ■ q{p2 ■ P3P1 ■Pi+P2- PiPi ■P3-Pi- P2P3 ■ Pi) 



mr, 



(A4) 



for Ti — > VIV decays. The pi, p2, P3 and p4 are 4- momenta of a heavy meson (7i), a hght 
meson (P of V^), a lepton (/) and antilepton (f), respectively. The corresponding masses 
are mi, r?7,2, ""^3 and 777.4. The phase functions contain integration over Mandelstam variable 
"S13 = (Pi ~ P3Y- The limits of integration are defined in the standard way |^ . 
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Table I. List of neutrinoless LFV processes, the composite loop functions and the tree 
level functions contributing to them and the approximations (physics) needed for evaluation 
of amplitudes. /, P, V, Ti and B denote leptons, light pseudoscalar mesons, light vector 
mesons, heavy pseudocsalar mesons (containing c or 6 quark), and light baryons, respectively. 
In the first column the list of the neutrinoless LFV processes is given, with references only 
to the calculations made within extensions of SM with heavy neutrinos. The abbreviations 
cqf = conserved quark flavor, ncqf = nonconserved quark flavor and H = Higgs mediated 
process, serve to distinguish processes with seemingly similar particle content. In the second 
column, the Feynman diagrams contributing to any specific process are listed. For instance, 
l-q-hox corresponds to the box diagram with one lepton current and one quark current. In 
the third column the approximations and physics used for calculation of amplitudes are 
listed. Following abbreviations are used: HQET = heavy quark effective theory, CHPT = 
chiral perturbation theory, VMD = vector meson dominance, GTR = Goldberger-Treiman 
relation, 1 = lepton physics, q = quark physics. 



process 


diagrams 


approximations (physics) 


/-/'7111a 


7 




/i ^ e conversion 1^,7^,71] 


7, Z and l-q-hox 


1, q, nuclear 


M ^ M conversion 


l-hox 


1, atomic 


/- ^l'-l{l^ 11,12,18 


7, Z and l-hox 




T -^ IP° (cqf) 11,|18| 


7, Z and l-q-hox 


1, q, PCAC 


r -^ IP° (ncqf) 11 


l-q-hox 


1, q, PCAC 


r -^ IV° (cqf) TT 


7, Z and l-q-hox 


1, q, VMD 


r -^ IV^ (ncqf) 11 


l-q-hox 


1, q, VMD 


z^ii' 11J63 


7, Z and l-hox 




H -^ ir |es|§| 


H 




PO ^ e/. (cqf) 


7, Z and l-q-hox 


1, q, PCAC 


PO ^ e/i (ncqf) 10 


l-q-hox 


1, q, PCAC 


n'^ -^ //' 


l-q-hox 


1, q, PCAC 


r - /'-P1P2 (cqf) m 


all except l-hox 


1, q, CHPT, PCAC, VMD 


r- -> /'-P1P2 (ncqf) [|| 


l-q-hox and W^W~ 


1, q, CHPT, PCAC, VMD 


T- -^ /'-P1P2 (cqf,H) H 


H and W+W- 


1, q, CHPT, PCAC 


Pi -^ Pse/i ITO| 


l-q-hox 


1, q, VMD, CHPT 


n -^ PW 


l-q-hox 


1, q, VMD, CHPT, HQET 


n -^ vir 


l-q-hox 


1, q, VMD, CHPT, HQET 


Pi -^ B2efi 10| 


l-q-hox 


1, q, PCAC, GTR 
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Table Ila. Coefficients defining the meson content in axial-vector quark currents with 
denoted quark content and normalization given by Eq. (0). Two additional coefficients are 



J -a: -t c boxAs i j box,sd -i 

ditierent irom zero: a^o = 1 and a^o = 1- 



pO 


a|o 


box,uu 

apo 


box Ad 
apo 


apo 




-V2 

\/2cp Sp 

cp \/2sp 
V3 V3 


1 

cp 1 Sp 

a/6 ^ a/3 
Cp sp 
a/3 a/6 


1 

cp Sp 

a/6 a/3 

cp _|_ sp 
a/3 ^ a/6 




a/2cp Sp 
V3 a/3 
cp V2sp 
a/3 a/3 



Table lib. Coefficients defining the meson content in vector quark currents with denoted 
quark content and normalization given in Eq. (0). Two additional coefficients are different 



from zero: a 



box,ds 
K*0 



—1 and a 



box,sd 
K*0 



-1. 



v 


a^o 


box,uu 

ayo 


box Ad 

ayo 


6ox,ss 


/37 


7yo 


UJ 


V^CvC2W 1 Sy 

V3 '+ V3 
cv 1 \/2svC2w 
a/3+ a/3 


1 

cv sv 

Cv _|_ Sv 
a/3 + a/6 


1 

Cv _|_ Sv 
a/3 a/6 




V2cv' 1 Sv 

Cv 1 \/2sv 
a/3+ a/3 


2a/2cvs|, 

V3 

2V2svs|, 

a/3 


2V2cvs2^, 

V3 

2V2svs2j, 

a/3 
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Table Ilia. The comparison of experimental and theoretical upper bounds on LFV BRs. 
Experimental upper bounds for unmarked processes are taken from Ref. 
denoted by * are from Ref. |^. The newest value B^^{^~ -^ e~7) = 1.2 x 10 
in Ref. 



while those 
^^ is given 
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Process 


^exp 


r>UB 


#^- 


^e 7 


3.8 X 10-1^ 


8.08 X lO-^x^^ 


r — 


> e"7 


2.7 X 10-6 


3.38 X lO-^x^g 


r^ — 


*/""7 


3.0 X 10-6 


6.68 X 10-9a;2^ 


^^ - 


* e^e+e~ 


1.0 X 10-12 


6.41 X lO-^yl^ 


r^ — 


> e"e+e~ 


2.9 X 10-6 


2.69 X 10-6y2^ 


r^ — 
r^ — 
r^ — 




1.9 X 10-6 
1.8 X 10-6 
1.7 X 10-6 


4.48 X 10-^y2^ 
1.44 X 10-6y2^ 
3.71 X 10-^y2^ 


r^ — 
r^ — 
r^ — 




1.5 X 10-6 
1.5 X 10-6 
3.7 X 10-6 


1.32 X 10-9y2^^^ 
6.67 X lO-'Y,,,^ 
2.77 X 10-6y2^ 


r^ — 


>/^-7r0 


4.0 X 10-6 


5.40 X 10-^y^^ 


r^ — 


> e^?7 


8.2 X 10-6 


4.01 X lO-^y^^ 


T^ — 
T^ — 




9.6 X 10-6 
2.0 X 10-6 


7.81 X 10-8^/2^ 
2.70 X 10-6^/2^ 


T^ — 
T^ — 


> e 


6.3 X 10-6 
6.9 X 10-6 


5.27 X 10-^1/2^ 
2.30 X 10-6|/2^ 


T^ — 
T^ — 




7.0 X 10-6 
2.2 X 10-6 


4.46 X 10-^1/2^ 
2.67 X 10-6^/2^ 


T^ — 


> /X^TT + TT^ 


8.2 X 10-6 


5.19 X 10-7y2^ 


T^ — 


> e-K+K' 


6.0 X 10-6 


1.07 X 10-6^/2^ 


r^ — 


> fi-K+K- 


15 X 10-6 


2.07 X 10-^y2^ 


ttO- 


e~ ^^ 


1.72 X 10-s 


5.54 X 10-1^1/2^ 


Z ^ 




6 X 10-6 
1.7 X 10-6 


1.61 X 10-16^/2^ 
3.43 X 10-^y2^ 


#z- 


-^e^r± 


7.3 X 10-6 


8.08 X 10-6y2^ 


#z- 


-^/i^r± 


10 X 10-6 


1.59 X 10-6y2^ 


Process 


^exp 


T^UB 


*fi-Ti -^ e'Ti 


6.1 X 10-13 


1.01 X 10-5y2^ 


T~ — 


>e-Js:o 


1.3 X 10-3 


9.82 X 10-i6a;2^ 


T^ — 


>^-irO 


1.0 X 10-3 


1.93 X 10-16^2^ 


T^ — 


> e-ir*o 


5.1 X 10-6 


2.40 X 10-i5a;2^ 
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Table Illb. Theoretical upper bounds for some interesting BRs, for which experimental 
upper bounds have not been found. 
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FIGURES 

FIG. 1. The BRs and the upper bounds of the BRs for four leading processes of the 
groups of processes given in equations (p^), (|65D, ( ]66| ) and (^). Each of these four processes 
is shown in one of four panels. The BRs are evaluated for degenerate heavy neutrino masses, 
rriTv The model parameters are adjusted so that the maximal BR values are smaller than 
the present experimental upper bounds. It is assumed that the "angle" parameter x^r ~ 1, 
in accord with the Super-Kamiokande measurements. The full line represents the upper 
bound calculation keeping the parameters xw equal to one, and adjusting the s^ parameters: 
{sl^Y = 4.29 X 10-1°, (sl^f = 1.4 X 10-^ and {s^y = 3.3 x lO^^. The heavy long-dashed 
line represents the upper bound calculation keeping the s^ equal to the present experimental 
upper bound values (|), and adjusting the Xui parameters: x^e = Xre = 2.459 x 10"^, x,-^ = 1. 
The dotted line represents the calculation with the "realistic" Bin matrix elements, and with 
the same parameters as for the full line calculation. 



FIG. 2. The BRs and upper bounds for BRs for the same processes as in Fig 1., but now 
evaluated as a function of the ratio of two heavy neutrino masses, mN^/mNi- For all curves 
the first and third masses are taken to be degenerate, m^-^ = rriN^ = 4000 GeV, while the 
second mass assumes values within the interval 1 < 771^2/ f^tNi < 10^. The types of lines 
represent the same sets of parameters as in Fig. 1. In the first panel, representing the /i — >■ e 
conversion on Ti, additional curve is added, to show that one can always achieve theoretical 
values smaller than the present experimental bounds. The calculation for that curve was 
made within the upper bound procedure, for xui = 1, (s^")^ = (s^'')^ = 0.5 x 10^^° and 
(sl^Y = 0.033. 
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